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In this Letter we report on a search for long-Hved particles that decay into final states with two 
electrons or photons. Such long-lived particles arise in a variety of theoretical models, like hidden 
valleys and supersymmetry with gauge-mediated breaking. By precisely reconstructing the direction 
of the electromagnetic shower we are able to probe much longer lifetimes than previously explored. 
We see no evidence of the existence of such long-lived particles and interpret this search as a quasi 
model-independent limit on their production cross section, as well as a limit on a long-lived fourth 
generation quark. 
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PACS numbers: 13.85.Qk,13.85.Rm,14.65.-q,14.80.-j 



The standard model is surprisingly successful in de- 
scribing phenomena observed at accelerators. One would 
expect, given its numerous theoretical shortcomings and 
the proliferation of searches for deviations from it, that a 
more general underlying theory would have been already 
revealed. It is therefore a possibility that the discov- 
ery of new physics eludes us because the new physics 
looks different from popular standard model extensions 
like minimal supersymmetry (SUSY). 

In this Letter we search for pairs of electromagnetic 
(EM) showers from electrons or photons that originate 
from the same point in space, away from the pp interac- 
tion point. Such events can be a signature of a long-lived 
b' quark decaying into a Z boson and a jet ^1]. In mod- 
els with gauge-mediated SUSY breaking [3| a long-lived 
neutralino with large higgsino component can decay into 
a Z boson and a gravitino. In the hidden valley mod- 
els 0, t'-mesons can decay into electron pairs. In all of 
the above examples, a significant imbalance in transverse 
energy can be present due to Z boson or v hadron decays 
into neutrinos or lightest supersymmetric particles (LSP) 
that remain undetected. 

A search for such long-lived particles at hadron col- 
liders was performed by CDF based on the recon- 
struction of lepton tracks from a secondary vertex. The 
sensitivity to large lifetimes in that search is limited by 
the difficulties in reconstructing tracks that originate far 
from the interaction point. In our analysis, we use the 
fine segmentation of the DO detector to reconstruct the 
directions of the EM showers and use that to reconstruct 
the common vertex. This method allows us to probe 
dramatically longer decay lengths, albeit at the price of 
lower sensitivity to short lifetimes. Since we do not re- 
quire the electron track to be reconstructed, our search 
results are also applicable for long-lived particle decaying 
into photons. 

The data in this analysis were recorded with the DO 
detector [HI, which comprises an inner tracker, liquid- 
argon/uranium calorimeters, and a muon spectrometer. 
The inner tracker is located in a 2 T superconduct- 
ing solenoidal magnet and consists of silicon microstrip 
and scintillating-fiber trackers. It provides measurements 
of charged particle tracks up to pseudorapidity Q of 
?7| w 3.0. The calorimeter system consists of a central 
section (CC) covering |?7| < 1.2 and two endcap calorime- 
ters extending the coverage to |?7| ~ 4, all housed in sep- 
arate cryostats 0. The electromagnetic section of the 
calorimeter has four longitudinal layers and transverse 
segmentation of 0.1 x 0.1 in 77 — space (where (p is 
the azimuthal angle), except in the third layer, where 
it is 0.05 X 0.05. The central preshower (CPS) system 
is located between the solenoid and the CC calorime- 
ter cryostat, covers |?7| < 1.2, and provides measurement 



of EM shower position with a precision of about 1 mm. 
The data for this study were collected between 2002 and 
summer 2006 using single EM triggers . The integrated 
luminosity Q of the sample is 1100 ± 70 pb~^. 

We select events with two EM clusters reconstructed 
in the central calorimeter with transverse momentum 
Pt > 20 GeV and \r]\ < 1.1, with the shower shape con- 
sistent with that expected of a photon. EM clusters are 
required to be isolated in the calorimeter and tracker . 
Both EM clusters are required to have a matched CPS 
cluster. Only CPS clusters within a fixed rj — (j) window 
are considered for matching, limiting the electron or pho- 
ton distance of closest approach {DCA) to the beam line 
at approximately 16 cm. Jets are reconstructed using the 
iterative midpoint cone algorithm [lo| with a cone size of 
0.5. The missing transverse energy is determined from 
the energy deposited in the calorimeter for jr;! < 4 and is 
corrected for the EM and jet energy scales. 

The DO EM pointing algorithm fits five shower posi- 
tion measurements (one in the CPS and four in the four 
EM layers of the central calorimeter) to a straight line 
which is assumed to be the EM object direction. The 
electron trajectory for energies above 20 GeV, which are 
of interest to this analysis, is very close to a straight 
line, which is defined by the energy-weighted EM cluster 
position [x^'^^.y^^^) and the DCA. The DCA recon- 
struction accuracy is about 2 cm. The common vertex 
position in the xy plane for two EM objects is the inter- 
section of the two lines associated to them and is given by 
a solution of the system of two linear equations (see Fig. 
[l]for definitions of the trajectory and quantities below): 

/ -Ayi Axi \ / x \ / • Axi - xf • Ayi \ 

1^ -Ay2 Ax2 j U / V y2^'' ■ - x§^^ -Ay^J- 

The determinant of this system, D, is proportional to 
the sine of the opening angle 6*12 between the EM ob- 
jects. The vertex transverse position resolution is in- 
versely proportional to the determinant. Therefore, in 
the following we consider events with \D\ > 4000 cm^ , 
which roughly corresponds to sin0i2 > 0.5, and use the 
variable Rs = ii/x^ + y"^ ■ (D/1000 cm^), which, while 
related to the reconstructed vertex radius, also takes into 
account its uncertainty. The sign of Rs is given by the 
sign of the scalar product of the px of the pair of EM 
objects with the vector pointing from the origin to the 
vertex location of the two EM particles. To reduce the 
background we further require that at least one of the 
two EM objects has DCA > 2 cm. 

For vertices that originate from real particle decays, 
Rs is positive, while its distribution for prompt electron 
or photon pairs is symmetrical around zero. The latter 
assumption was extensively checked with Monte Carlo 
(MC) simulation, a Z ^ e"*"e~ data sample (both elec- 
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TABLE I: Observed number of events {Rs > cm) and esti- 
mated background {Rs < cm) for different selections. 



FIG. 1: Definition of the reconstructed EM particle trajec- 
tory. In the DO coordinate system the equation of the trajec- 
tory is given by Ax ■ {y ~ y^"^^) = Ay ■ {x ~ x^^^). The 
distance from the beam line to the EM shower maximum 
^ (a;CALp _|_ (^yCAL'j2 jg typically around 90 cm. 



irons in the Z — > e+e~ sample were required to have 
reconstructed tracks originating from the primary ver- 
tex), and a control sample of multi-jet events which has 
been selected exactly as the signal events except with 
an inverted tracker isolation requirement. Therefore, we 
estimate the background for positive values of Rs by mir- 
roring the negative part of the distribution. 

The invariant mass M of the two EM objects is cor- 
rected for the reconstructed vertex position, and the data 
are divided into three bins: 20 < Af < 40, 40 < M < 75, 
and Af > 75 GeV. The last bin is used for searches for 
the fourth generation 6'. The corresponding observed Rs 
distribution is shown in Fig. [2l All mass bins are used 
for a quasi model-independent search for long-lived par- 
ticles. We also examine events with > 30 GeV and 
Af > 20 GeV. No excess of events with positive Rs val- 
ues is present in data (see Table so we proceed to set 
limits on new physics. 

_We use PYTHIA 6.319 [111] to generate events pp 
b'b' ZbZb -> e+e- -I- X. PYTHIA calculates pro- 
duction cross sections varying from 79.4 to 3.6 pb as 
the b' mass changes from 100 to 190 GeV. The events 
are then processed through the GEANT-based MG 
simulation, electronics and trigger simulation, and are 
reconstructed with the same reconstruction program as 
collision data. The expected Rs distribution for a typical 
signal point is shown in Fig. [21 We use the efficiencies 
and acceptances obtained using this signal MG for the 
model-independent search as well. The significant jet 
activity in these events gives a conservative estimate of 
the efficiency for SUSY scenarios and should be adequate 
for hidden valley models [l3|- In order to study differ- 
ent masses of hypothetical resonances in addition to the 
samples above we also generated samples of b' — > vb for 
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V masses of 30 and 50 GeV. We find that the efficiency 
and acceptance for the MC events have no significant de- 
pendence on the masses of the 6' and v. We set the b' 
mass to 150 GeV and vary its lifetime cr between 2 and 
7000 mm. 

In Fig. [3] we display the limits on the production cross 
section of a long-lived particle times its branching frac- 
tion to decay into a pair of electrons. Limits were ob- 
tained from the Rs distribution using the modified fre- 
qucntist approach ^l3| as implemented in [l5l |. This 
method is based on a log-likelihood ratio [LLR) test 
statistic, and involves the calculation of confidence lev- 
els for the signal plus background and background-only 
(null) hypotheses (denoted by CLg+t and CLb, respec- 
tively) by integrating the LLR distributions resulting 
from simulated pseudo-experiments. The upper limit on 
the cross section at the 95% C.L. is defined as the cross 
section value for which the ratio CLg = CLg+b/CLi, — 
0.05. The systematic uncertainties were taken to be flat 
as a function of Rs- They include the uncertainty in 
electron or photon identification and triggering (15%), 
uncertainty on Monte Carlo simulation (5%), and un- 
certainty on luminosity (6.1%). At the cr value of 100 
mm we exclude at the 95% C.L. the production cross 
section times branching fraction of long-lived particles 
that decay into a pair of electrons or photons above 1.9 
pb, 10.2 pb, 7.1 pb, and 4.4 pb for > 30 GeV and 
Af > 20 GeV, 20 < Af^ < 40 GeV, 40 < Af < 75 GeV, 
and Af > 75 GeV, respectively (see Fig. [3]). 

Intersecting the cross section upper limits shown in 
Fig. 3d with the theoretical cross section of the produc- 
tion of the fourth generation b' quark IjJ we compute 
limits on its lifetime as a function of its mass assuming it 
decays only into Zb. The limits are presented in Fig. 31 
together with the exclusion region from the track-based 
CDF search H] . The two search methods are complemen- 
tary to each other. 

To summarize, we have performed a search for long- 
lived particles decaying into electron or photon pairs us- 
ing a new method that allowed us to explore previously 
unreachable portions of the parameter space. We find 
no evidence for such particles and present the results as 
model-independent limits on their production cross sec- 
tion and interpret them in the framework of a model with 
a long-lived b' quark [1]. 

We would like to thank Matt Strassler for many fruitful 
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FIG. 2: Observed Rs distribution for di-EM pairs with mass 
greater then 75 GeV (black points), expected distribution 
from prompt sources with its uncertainty (shaded rectangles) 
and the expected distribution in presence of b' quark with 
mass of 130 GeV and lifetime cr — 300 mm (solid line). 
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FIG. 3: Expected (white triangles) and observed (black 
squares) 95% C.L. upper limits on the cross section of a long- 
lived particle times the branching fraction of its decay to ei- 
ther a pair of electrons or photons for (a) ]^rp > 30 GeV and 
M > 20 GeV, (b) 20 < M < 40 GeV, (c) 40 < M < 75 GeV, 
and (d) M > 75 GeV. All observed upper limits are within one 
standard deviation (shaded band) from the expected limits. 
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FIG. 4: 95% C.L. exclusion region of b' lifetime (cr) vs. mass 
for CDF Run I and current DO result. 
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